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1. Introduction 

The nature of the dark matter (DM) in the Universe remains one of the most challenging 
problems in cosmology. Numerous candidates for DM have been proposed in the literature. 
One of the most popular candidates in the context of supersymmetric theories with i?-parity 
conservation is the lightest supersymmetric particle (LSP). Such a candidate has usually 
been studied in the framework of supersymmetric extensions of the Standard Model, such 
as the Minimal Supersymmetric Standard Model (MSSM) |l[ or its constrained version 
(CMSSM) P]. In these models, the LSP is often the lightest neutralino, i.e. a linear combi- 
nation of the fermionic superpartners of neutral gauge and Higgs bosons. The interactions 
of the neutralino are weak, and its number density at decoupling is therefore often of the 
required order of magnitude, which makes it an excellent candidate for the WIMP (Weakly 
Interacting Massive Particle) ||. 

The case of the neutralino WIMP is very appealing and puts strong bounds on SUSY 
masses and parameters but it is certainly not a unique one. One should not forget 

that very different possibilities can arise in other well-motivated extensions of the Standard 
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Model. For example, in the framework of supergravity there arises the gravitino, the SUSY 
partner of the graviton. The gravitino has long been known to be a potential candidate for 
DM 0], A number of more speculative possibilities, stemming from, for example, string 
theories, have also been proposed. 

The aim of this paper is to investigate another well-motivated scenario, in which the 
LSP is the axino. Such a particle is present in SUSY models that incorporate the Peccei- 
Quinn (PQ) solution to the strong CP problem. The PQ global U(l) symmetry || was 
introduced more than 20 years ago to solve the strong CP problem, and still remains the 
most elegant solution, especially since lattice results seem to disfavor the possibility of a 
massless up-quark Q. The axino does not belong to the usual MSSM spectrum, but, as 
a fermionic superpartner of the axion ]?], ||, is much more weakly interacting. Similarly 
to the case of the axion, the couplings of the axino to "ordinary" matter (below simply 
referred to as "matter") are suppressed by the Peccei-Quinn scale f a ~ 10 11 GeV, thus the 
axino remains invisible to present experimental searches. 

A stable relic axino has been shown to be an attractive candidate for either warm ]9|, 
1C, 11, [T^] or, preferably, cold dark matter (CDM) 13]. An underlying assumption 



is that a primordial axino population was washed away by an earlier inflationary period. 
Axinos can then again be produced in large numbers through several possible mechanisms. 
One efficient way is through thermal production (TP); this proceeds via thermal scatterings 
involving, e.g., gluinos (and their decays) in the primordial plasma |To|, [l2|], in close analogy 
with gravitino production ||]. Another one is non-thermal production (NTP): an example 
of this being an out-of-equilibrium decay of the next-to-lightest particle (NLSP) |], A 
number of papers have investigated cosmological properties of axinos as warm |9|, O, 14] 
or cold relics (ITol [iU . 



In a recent extensive analysis [12], the role of the effective dimension-5 gluino-gluon- 
axino interaction was examined in the production of stable axinos as relics in the Universe. 
At high reheat temperatures Tr (much larger than squark/gluino masses but still Tr <C f a 
in order not to restore the PQ symmetry), the thermal production of axinos via scattering 
processes was shown to be dominant. This allows one to place strong bounds on the reheat 
temperature as a function of the axino mass by imposing the constraint £lah? < 1. At 
lower Tr, gluino decays in the thermal plasma and, at Tr <C mj, non-thermal production 
of axinos were shown to become more important. 

The scope of this paper is to supplement the analysis in |l3|, [l2| by investigating the 
thermal production of axinos involving the dimension-4 squark-quark-axino interaction. 
The effective vertex for this interaction arises at a one loop order in an effective theory but 
we find it to be enhanced by the gluino mass. For this reason, we argue that squark decays 
cannot be neglected in a computation of the axino abundance, especially in the case of a 
heavy gluino. 

The paper is organized as follows. In Sec. || we briefly review axion models. In Sec. || 
we introduce the effective squark-quark-axino vertex which is then computed in Sec. |||. 
The formalism for the thermal production of axinos is reviewed in Sec. || and applied to 
two- and three-body squark decays in Sec. [6| and to scattering of quarks and squarks in 
Sec. ^. We present our numerical analysis in Sec. |8| and conclude in Sec. ^. 
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2. Axion models 



In the PQ scenario, a complex scalar field is introduced, which breaks the global U{V}pq 
at a high scale and whose Goldstone boson component, the axion, plays the role of a 
dynamical Oqcd relaxing at the origin after the QCD phase transition. In fact, after the 
chiral symmetry breaking, the axion acquires a tiny mass from instanton effects [17|. 



Two classes of models implementing the PQ mechanism are usually studied in the 
literature: the KSVZ-type (hadronic) axion models (KSVZ) 111] and the DFSZ-type mod- 



els [19]. In the first case, the SM particles are not charged under U(1)pq, while in the 
second they are. For this reason, different axion interactions are present in the two ap- 
proaches. 

In supersymmetric extensions of axion models the axion field is promoted to a chiral 



multiplet [20 1 which contains not only the pseudoscalar axion, but also a scalar, the saxion, 
and their fermionic superpartner, the axino. An interesting feature of such supersymmetric 
models is that the PQ global symmetry is enlarged to a complex U(l), which ensures that 
the whole axion multiplet is massless and degenerate pl| , as long as supersymmetry and 
the chiral symmetry are unbroken. If supersymmetry is broken at low energies, as is 
required to explain the large hierarchy between the electroweak and the Planck scale, we 
can assume that this degeneracy is not strongly lifted, and consider a model in which the 
axion multiplet is the only additional low-energy degree of freedom compared to those 
already contained in the MSSM. 

The breaking of supersymmetry reduces the complex U(l) to the usual U(1)pq and 
generates a soft mass for the saxion [22], while the axion remains protected by the surviving 
U{\)pq. For the axino different scenarios can arise, strongly dependent on the model and 
the particle content [^, 23]. In general, the axino, being a component of a chiral multiplet, 
does not acquire a tree-level Majorana soft mass term. However, a Dirac mass term may 
arise which will lead to mixing with the neutral heavy states in the PQ sector and the 
MSSM neutral higgsinos, giving rise to an enlarged neutralino mass matrix at low energy. 

On the other hand, if the mixing of the neutralino-axino is very tiny, or absent, an axino 
nearly massless mass eigenstate survives at the tree level (see [21, 23, 12]). Particularly 



interesting is the case in which a Majorana mass for the axino is generated at one loop [24, 



25], falling naturally in the range of a few tens of GeV, often below the present bounds on 
the neutralino. This bound obviously does not apply to an axino state due to its suppressed 
interactions with SM particles. 

Phenomenologically, an appealing scenario arises: the LSP axino remains unobserved 
due to the very tiny couplings of the axion multiplet with matter, and does not spoil the 
precision measurements of the SM, but still plays a key role in the cosmological evolution 
of the Universe. 

3. The interactions of the KSVZ axino 

In this paper we consider mainly an axino of the KSVZ type; our final results are expected 
to be similar to the case of the DFSZ axino, and we will comment on the differences below. 
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In the KSVZ, or hadronic, axion models, the axion multiplet (and therefore also the 
axino) does not couple directly with ordinary matter, instead its interactions with matter 
are suppressed not only by the PQ scale, but also by loop factors. 



A simple example for the superpotential for a KSVZ axion multiplet is [26| 

\$QQ + S(f% 



W 



(3.1) 



where Q,Q are PQ-charged, heavy colored triplet states, $ and <£> are PQ-charged SM 
singlets, while S is a singlet with respect to both the SM and PQ symmetries. 

The axion multiplet is in this simple case contained in $ and writing their scalar 
and fermionic part in terms of the axion multiplet A = (A, a) and another complex scalar 
and fermion pair (B,ip), we obtain respectively 



/ a exp 
/ a exp 



(B + A)/(V2f a )\ , (i> + a)/y/2j 
(B - A)/(V2/«)] , (V - a)/V2) • 



(3.2) 



From the superpotential (|3.1|) , we see that in the supersymmetric vacuum with = 
f%, only the axion multiplet remains massless and all the other states acquire a mass of 
the order of f a . Note also that no direct coupling between the axion and the MSSM fields 
is present in ( |3.1| ). However, axino couplings to MSSM fields do arise in the low energy 
effective action and will be derived below. 





Figure 1: The two Feynman diagrams contributing to the gluino-gluon- axino interaction. The 
extra quark Q and squark Q states are all very heavy, Trig q ~ f a , which allows one to contract the 
loops to an effective "tree-level" gluino-gluon-axino vertex. 

The interaction of the axino with gluons and gluinos proceeds via the anomaly one-loop 
diagrams shown in Fig. |l|. After integrating out the heavy KSVZ quarks Q and squarks Q 
(rrig g ~ f a ), we are left with a single effective gluino-gluon-axino dimension-5 interaction 
term in the Lagrangian 



a l5 a^g b G b 



[IV 



(3.3) 
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where b = 1, . . . , 8, g is the gluino and G is the strength of the gluon field, while N is the 
number of flavors of quarks with the Peccei-Quinn charge, and is equal to 1 (6) in KSVZ 
(DFSZ) models. For the remainder of this paper, we understand f a to mean f a /N. 

The coupling fl3.3| ) can also be obtained by supersymmetrizing the axion-gluon-gluon 
interaction, which corresponds to adding the following Wess-Zumino term to the MSSM 
superpotential: 

Wwz = 1^Y A Tr [W a W a ] , (3.4) 

where A is the axion multiplet, W a the vector multiplet containing the gluon, and the trace 
sums over color indices. An interesting characteristic is that the coupling above is only de- 
termined by the QCD anomaly of the heavy states and is not subject to renormalization [27] 
or model dependence. 

In an analogous way, Wess-Zumino terms can arise also for the other gauge interactions, 
depending on the charge of the heavy quark multiplet; at high energy, when all leptons can 
be considered massless, such interactions can be rotated into the U(1)y direction |12| and 
we are left to consider only the term 

W wz = a J ( l aY J AB a B a , (3.5) 



where B a is the hypercharge vector multiplet and C a YY is a model-dependent factor [28], 
which vanishes if the heavy KSVZ quarks are electrically neutral. If they have electric 
charge e<g = —1/3, +2/3, then C a yy = 2/3, 8/3. This superpotential term leads to the ad- 
dition of the following effective dimension-5 interaction term to the low-energy Lagrangian 

O-yCoYY- 

8vr/ a 

Cosmological implications of the effective axino-gauge boson-gaugino operators (|3.3|) and 



C aBB = LT hscr^BB^ (3.6) 



( |3.6| ) have been extensively studied in [12]. 

However, in addition there arises the effective dimension-4 coupling of the axino to 
quarks and squarks when we integrate out the heavy states: 

£aqq = 9^ff <Lj QjP R/Lj 5 ®- (3.7) 

We find that this effective coupling, which arises at two-loop level in KSVZ models (and 
therefore at a one-loop level in the effective theory valid much below f a ), can play an 
important role in the thermal production of axinos. 

Instead of undertaking the complicated task of a full computation of g e ^ in the renor- 
malizable theory, which would give a result strongly dependent on the unknown parameters 
of the axion superpotential anyway, our strategy is to resort to the expedient of estimating 
the coupling via an effective theory, given by the MSSM plus the interaction in Eq. (p.3j). 
That is to say, we compute the two one-loop diagrams in Fig. and use our result to define 
an effective squark-quark-axino vertex in the low energy Lagrangian. We are then able to 
study the effects of the coupling on the axino abundance. 

Let us now mention the situation for the DFSZ models. In that case, the axino can 
mix with all the other neutralinos through the coupling which generates the Higgs [x term, 
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Figure 2: The Feynman diagrams contributing to the squark-quark-axino interaction. The thick 
dot denotes the effective gluino-gluon-axino vertex of Fig. [I]. 

and due to this mixing the coupling in Eq. ( |3.7|) arises at tree level. Since the mixing is 
strongly suppressed by at least the ratio of the electroweak to the PQ scale, vew/ fa, anci 
since we consider the case when an axino-like light state remains in the spectrum, we do 

T I R 

not expect the full DFSZ result for g e jj- to be much larger than our result for the KSVZ 
case. 

In general, the effective coupling to a quark and a squark could be enhanced for 
the DFSZ axino, but is not expected to be suppressed, unless there are very unlikely 
cancellations between diagrams. Therefore, the upper bounds on the reheat temperature 
we obtain for the KSVZ case are also conservative bounds for the DFSZ model. 



4. The effective squark-quark-axino vertex 

We now proceed to compute the effective dimension-4 squark-quark-axino interaction de- 
fined by the two diagrams in Fig. |2[ The corresponding amplitudes for an axino a with 
momentum k\ , a quark q with momentum and a squark q L l R with momentum k\ + hi , 
are given by: 

A fR = ^jf'VifcOfei) (4.1) 

<J la 

A' R = T^f u 9 (W 2 L/ V^i) (4.2) 

where summation over color has been performed, the upper/lower signs refer to left- and 
right-handed (s)quarks respectively, and I L < R, s come from loop integrations. We have 
suppressed the flavor index since the amplitudes are diagonal in flavor space. 

J / R 

The integral of I x is logarithmically divergent; in the full renormalizable KSVZ 
model (with broken supersymmetry), such a divergence would be taken care of by a coun- 
terterm and reabsorbed in the parameters of the Lagrangian. Rather than perform the full 
renormalization, we instead keep the dependence on the UV cut-off f a . This procedure is 
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justified by the fact that we are working in the context of a non-renormalizable effective 
theory, valid only below f a . At higher scales, all the heavy degrees of freedom become 
dynamical. With this method, we are absorbing all our ignorance about the high energy 
theory into one single parameter, the "effective" PQ scale f a . This allows us to perform a 
model-independent analysis of the implications of the effective coupling. 

We have performed the loop integration using the usual methods, based on the Feyn- 
man parameterization. To lowest order in l// a , for m qj — ► and assuming that the axino 
mass is much smaller than the gluino mass rrig, we obtain 



T L/R 



T L/R 



Sirrig 
Sirrig 



log 



fa 



111; 



111'- 



+ 



in 



2k x ■ k 2 



log 1 



2ki • k 2 



(4.3) 



1 TT Z 
1 »_ 

3 18 2fci • k 2 



1 



3 2ki • k 2 



-Lb 1 



2k x ■ k 2 



mf 



+ 



1 
6 



1 + 



m- 



2ki ■ k 2 



log 



2ki ■ k 2 



in- 



(4.4) 



Due to the chiral structure of the loop, the result is proportional to the mass of the internal 
fermions. 

Using these expressions and the definition of g^lf in (3/7), we arrive at 



L/R 
9eff 



=R- 



AV2 



- =F 



fa 
171; 



T L/R 



2I^ /R ) 



v^vr 2 fa 



log 



fa 



(4.5) 
(4.6) 



Note that the dominant contribution comes from the logarithmically divergent part. 
However, as in all supersymmetric theories the dependence on the cut-off is weak. Also, 
the limit of large f a is perfectly well-defined and the coupling vanishes in the decoupling 
limit, when f a — > +oo. Another important point is that the coupling depends on a single 
unknown parameter, the ratio m g /f a : therefore in the following we keep f a = 10 11 GeV 
fixed, and vary the value of the gluino mass to explore the parameter space. 

Secondly, the effective vertex is proportional to the gluino mass and tends to zero 
in the limit of exact supersymmetry and massless quarks; this can be easily understood 
by considering that for massless quarks the whole axion multiplet decouples from the 
low energy Lagrangian [H]. When the chiral symmetry is broken, a term proportional to 
the quark mass survives, giving an axino-quark-squark coupling of the same type as the 
standard axion-quark-quark coupling. 

However, in the case of broken supersymmetry a new contribution to the coupling 
arises from the gluino mass and dominates for all light quarks. For the top quark, the term 
proportional to nit could be important for small gluino masses. We neglect this term here 
since it would not change our results substantially. 

It is also worth mentioning that the full expression above contains, in the limit of 
negligible final state masses, terms proportional to log(l— m|/m|). Therefore g^jf contains 
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an imaginary component if mq > rrig. This imaginary piece corresponds to the "cut" loop 
diagram, in which the intermediate gluino and quark are real. Since this contribution is 
already included when considering the squark decay into a gluino and a quark, we keep 

T I Ft 

only the real part of g e ^ to avoid double counting. 
5. Thermal production of axinos 

For reheat temperatures below f a , the axinos are not in thermal equilibrium with the MSSM 
particles. Assuming that any primordial axino population was completely diluted away 
during an inflationary stage, the present axino density has been generated after the end of 
inflation, since the epoch of reheating/preheating. Two different production mechanisms 
are possible: thermal production, TP, in which axinos are produced via thermal processes 
in the primordial plasma; and non-thermal production, NTP, which is independent of the 
presence of the thermal bath. An example of the latter is the decay of the neutralino NLSP 
after it freezes out, as considered in [13]. In this case the final number density is related 



only to the density and lifetime of the initial neutralino population. 

In the case of TP, on the other hand, the axino number density n„ depends on the 
species and processes present in the plasma, and can be obtained by integrating the Boltz- 
mann equation with the appropriate collision integral. Fortunately, since the number den- 
sity of axinos is well below the equilibrium one for Tr <C f a , we can neglect the inverse 
processes, which are suppressed by n„, and we have: 

+ 3Hn & = y^(r(i -» a H ))ui + y~]((r(i + j a H )v ie i)mnj, (5.1) 

i i,j 

where the Hubble parameter, H, can be given as a function of the temperature by H(T) = 

1 /2 

[{^g*) I (90M|)] 1 T 2 in the radiation-dominated era, with counting the number of 
effective massless degrees of freedom and Mp denoting the reduced Planck mass, Mp = 
2.4 x 10 18 GeV. 

The first term on the r.h.s. takes into account all the decay channels of the ith particle 
with one axino in the final state, and (T(i — > a + • • •)) is the thermally averaged decay 
width for the process. The second term corresponds to axino production via 2 to 2 particle 
scatterings, with o~(i + j — * a + • • •) denoting the cross-section for particles i,j scattering 
into final states involving axinos and v re \ the relative velocity. Averaging over initial spins 
and summing over final spins is understood. 

In order to solve Eq. ( |5.1| ), it is convenient to separate the two contributions and 
introduce the axino TP yield as p 



,,,TP 



/TP n « 



a 



<£ Y *« + Y I Y% A , (5.2) 



where s = (2tt 2 /45)g s *T 3 is the entropy density, and normally g Sif = in the early 
Universe. 
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Moreover, by changing variables from the cosmic time t to the temperature T, we can 
write the two solutions of the Boltzmann differential equation easily in integral form: 



Y**(T ) = QiT ^^^ (5.3) 

where we have considered the evolution from the reheating temperature after inflation, Tr, 

hi 



down to the present temperature To. Full expressions for Yf™ and K s f * can be found 



m 



m 



In the following we consider the axino yield generated by the effective coupling in 
Eq. (|3.7| ). We compute the contribution due to squark decays and estimate the additional 
contribution to quark scattering. The axino number density arising from gaugino decays 
and scattering via the dimension-5 gluino-gluon-axino operator in ( |3.3|) was computed 



m 



6. Squark decay 

J / Ft 

The squark can decay into axinos in two different ways: via g e jj- into a two-body final state 
of an axino and a quark, or via an exchange of an intermediate gluino/neutralino and the 
anomalous couplings into a three-body final state of an axino, a quark and a gluon/photon. 
As we will see, the first one will play an important role while the second will always remain 
sub dominant. 



6.1 Two-body squark decay 

The process involves the diagrams shown in Fig. ||. The width for this decay of both left- 
and right-handed squarks is given by 



-\2body 



(9eff) 2 



16vr 

,4 



in,: 



a 



-m- n 



(2vr) 5 
~ (2.6 x 1(T 6 sec) 



^iog(A 

fa \m-g 



(6.1) 



ma 



500 GeV/ LVlTeV 



m g 



1 - 0.05 log 



m- n 



ITeV 



where in the last line we have taken the default value of f a = 10 11 GeV. This process is 
suppressed by a large power of the strong coupling and by loop factors, but it can still be 
sizable for large gluino mass, or equivalently for smaller f a . 



6.2 Three-body decay 

At an effective tree-level, the squark will also decay into an axino, a quark and a gluon, as 
illustrated in Fig. 



-9- 



Ql/r 




a 



Figure 3: Three-body decay of a squark into an axino, a quark and a gluon. The thick dot denotes 
the effective gluino-gluon-axino vertex of Fig. 1. 



The decay rate for this process in the case of large gluino mass rrig > rriq and negligible 
final state masses is: 

' " 1 (6.2) 



■^3body 



3(4vr)4 fi 



F -4 



in 



where we have summed over the color indices of the final state gluon and quark and where 



F(x) 



1 13x 7x 



6 



I + ~ + I (l-x)(3-7x)log 



{l-xf 



x- 



(6.3) 



In the limit m- g » m^, (|6.2|) reduces to 



-ySbody 



o; 



18(4^)4 m yi 



1 + \™§, 



(6.4) 



Note again that, in our effective theory, the 3-body decay width ( |6^ ) is proportional to 
the cube of the mass of the decaying particle. However, in the limit when one can neglect 
momentum flow in the intermediate states, one recovers in ( |6.4[ ) the usual fifth power of 
mass for a dimension-4 interaction. 
Comparing (|tT4]) with gives 



r 



3body 



TT 



r 2bod y i44 as 



in,. 



in; 



log" 2 (fa/mg) « 1 



(6.5) 



The 3-body decay channel is therefore subdominant for squarks lighter than gluinos. For 
example, if niq = 500 GeV, m q = 1 TeV and f a = 10 11 GeV, we have 



r 



3body 



(4.1 x 10" 2 sec) 1 



(6.6) 



which is suppressed by about four orders of magnitude relative to the two-body decay (|6.1|). 
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An analogous decay channel exists with an intermediate neutralino replacing the gluino 



and a final photon rather than a gluon. Taking the neutralino to be a bino [30], the decay 
of each color of squark into a quark, an axino and a photon is given by: 



a: 



( CaYyZ^Y- 



L/R 



(47r) 4 \ cos 9 W 



(6.7) 



L/R 



1/36,4/9,1/9 for a left-handed squark, a right-handed up-type squark 



where Yt 
q 

and a right-handed down-type squark, respectively. For m,g 
f a = 10 11 GeV, we have 



rf^aOsec)- 1 Cl YY Z\ x Yl 



L/R 



500 GeV, m„ = 1 TeV and 



(6.8) 



This is an even smaller contribution than that from the gluino exchange diagram above. 
Therefore we conclude that three-body decays are not expected to contribute substantially 
to the axino yield. 

Note that we have considered the decays only for the case of intermediate particle heav- 
ier than the initial squark. For an intermediate particle lighter than the initial squark, the 
dominant contribution to the axino production comes from the resonance. This contribu- 
tion can be factorized into the squark decay to a real neutralino/gluino, and its subsequent 
decay to an axino and a photon/gluon. This process is already included in the Boltzmann 
equation in the gaugino/neutralino decay term. 

6.3 Yield from the decay of squarks 

In general, the thermal average of a decay described by a width T is given by 



2vr 2 



m./T 



e x =F 1 



(6.9) 



where the — (+) is for decaying boson (fermion). By inserting this expression into (|5.3j), 
one obtains 



Y 



dec 



gMpT 
16ir 2 m 2 j tR 



dt- 



1 



e* =Fl 



7T 

2 " tan " 



+t R (t 2 - 2t\ 



t R 



(6.10) 



where t R = m t /T R and g = 135 v / T0/(27r 3 5 * /2 ), for g* = 915/4 for the MSSM in the early 
Universe. 

In the case T R > rrii, the integration gives the simple result 



ydcc ^ 3C(5)ffM P r 
1 \ixm 2 A 



3.11) 



- 11 - 




Q g 
Figure 4: Example of dimension-5 contribution to the process qq — > qa (Process I). 



To obtain the axino yield from squark decay, we substitute fl6.1| ) for the decay rate 
above, and obtain the following expression, valid for Tr > rriq 



Y; 



dec 



6C(5)c4 - fM P 
(2vr) 6 ' 



3 ( — 9 - 

rrin 



rrin 



Ydec ^ 



fa 



log 



fa 



10" 



(6.12) 
(6.13) 



The numerical value corresponds to m g = rriq = 1 TeV and f a = 10 11 GeV. It is clear that 
squark decay is typically as important as gluino decay in the thermal production of axinos, 
and can be substantially more important if the squarks are lighter than the gluinos. For 
Tr less than the mass of the decaying sparticle, both yields are Boltzmann suppressed and 
negligible. 

7. Scattering of quarks and squarks 

The effective dimension-4 squark-quark-axino coupling (|4.5| ) also gives new contributions 



to some of the scattering axino production processes previously considered in [12] which 
involved the effective gluino-gluon-axino vertex (|3.3| ). However, we shall see that such 
additional contributions will be always subdominant. 

As an example and for a conservative estimate, let us consider the scattering of two 
quarks into a gluino and axino. (This process was called process I in [J12|j.) This process 
does not contain sparticles in the initial state and is therefore never Boltzmann-suppressed 
by the number density of the initial particles. 1 Therefore process I is expected to be one 
of the dominant contributions to the axino yield from dimension-4 interactions at low 
temperatures. 

As we will show below, the dimension-5 diagrams that contribute to the scattering 
process, illustrated in Fig. |], are of order g^ |l2|| . The two dimension-4 diagrams that 
contribute are of order g^. They are shown in Fig. where the shaded squares denote the 
effective 1-loop squark-quark-axino vertices given in Fig. ^. Although these two diagrams 
are suppressed by g 2 s relative to the dimension-5 contributions, the enhancement through 
the large logarithm and the gluino mass present in g^ff partially compensates for it. 

1 Nevertheless, there is some suppression coming from final state gluino mass. 
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Generalizing the cross section due to the dimension-5 coupling computed in Ref. [12] 
by including a full dependence on the gluino mass, we find that 



07,5 



a' 



72vr2/ a 2 




(7.1) 



where the eight gluino species in the final state have been summed over, and the six flavors 
of quark in the initial state have been averaged over. 

The combined contribution from the two dimension-4 channels is 
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and is proportional to or s . The two squarks corresponding to each quark and the eight 
gluino species have been summed over in this expression, while we averaged over the initial 
flavor. 





g Q 
Figure 5: Dimension-4 contributions to Process I. 



Interference between the dimension-4 t and u channels gives a contribution proportional 
to ma, which is neglected here. 

There is also a contribution to the cross-section from interference between dimension-4 
and dimension-5 diagrams, which is of order a* . Giving the L squarks a common 
mass mi and the R squarks (qp) a common mass mjj, this 4/5 interference term can be 
written as 



°7,4/5 
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(7.3) 



Note that the contributions from the L and R squarks cancel in the limit = mp thanks 

J I Ft 

to the opposite sign of the coupling gl , . 

7.1 Yield from dimension-4 scattering processes 

We numerically computed the axino yield using an estimate of dimension-4 contribution 
to Process I derived above, and we found it to be negligible compared to the dimension-5 
part for any reheat temperature. 
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To illustrate this analytically, let us consider the following simplified example. In the 
limit vfiq = rrig = m, the contribution to the Process I cross-section from the dimension-4 
channel is 



07,4 
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(7.4) 



3s 2 

To calculate the yield we expand around m/Tp = 0. After summing over the four possible 
combinations of the initial spins, the six quark species and the three colors, and multiplying 
by a number density factor 3/4 for each initial fermion, we find 

Yft = 9 " s( ff ^ * ^ + O ( ) (7.5) 
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So for Tp greater than the masses of the external particles, the axino yields from these 
dimension-4 processes are approximately independent of Tp. For this reason, in analogy 
with decay processes, we have used a s (T = m) to take account of the running of the strong 
coupling. 

This is in contrast with what happens for the dimension-5 scattering processes. For 
example, the contribution of the diagram in Fig. |] to the yield is 

for high reheat temperatures, i.e. Y/5 oc Tr. This different behavior is simply due to 
dimensional reasons: the dimension-5 scatterings have cross-sections that are proportional 
to 1/ fa arid independent of s at large s, while the dimension-4 scatterings are characterized 
by a dimensionless coupling gHf oc rrig/ f a , and their cross-sections are suppressed by s _1 
at high energies. 

Comparing the two yields, we find that for Tr 3> m 

^r--/r(^) 5 (^) 3 (i) 

Therefore the contribution to the yield from the two dimension-4 diagrams is considerably 
less than the contribution from the dimension-5 diagram already considered for large re- 
heat temperature. As we will see, this is the case also for Tr ~ m, when the analytic 
formulas above are no more valid, since the dimension-4 yield reaches the asymptotic value 
in eq. ( |7.6| ) more slowly than the dimension-5 one. 

Below eq. ( [T.3|) we noted that the contributions to the cross-section from interference 
between dimension-4 and dimension-5 diagrams cancel in the limit mi = mp. However, in 
the case of unequal masses there is a remainder, and therefore a contribution to the yield 
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from 4/5 interference. We have calculated this contribution in two extreme cases; 

Y l%- \ Y tT ifrn L = Wm R (7.8) 

„J* Y *v* if mL = 0,i mR . (7.9) 

In either case, the magnitude of the effect of the dimension-4 diagrams does not become 
comparable with that of the dimension-5 diagrams. 

H 9 




'l 



Figure 6: An example of a 2-loop contribution to Process I. 

In principle, there are several other diagrams of order a\ that should be computed. For 
example, the box diagram in Fig. || will contribute to the total cross-section of Process I. 
Such diagrams are also 2-loop in a renormalizable theory (and therefore 1-loop in the 
effective theory) and as a result suppressed by the same loop factors, but contain no 
logarithmic divergence and so do not receive the enhancement factor contained in the 
diagrams considered above. Therefore it seems reasonable to neglect them too. 

We conclude that dimension-4 scattering diagrams never give a significant contribution 
to the thermal yield of axinos. The dominant contributions to Y~ cat are always those from 
dimension-5 diagrams already considered in |j~2]] . 

8. Results 

The relative contributions from the processes considered above to the axino abundance 
are summarized in Fig. |^. The solid line is the yield from dimension-5 scattering pro- 
cesses, while the dashed line is an estimate of dimension-4 scattering processes, computed 
by multiplying Yf± by seven, which is the number of channels receiving dimension-4 con- 
tributions. The figure illustrates that, for equal squark and gluino masses, the yield from 
dimension-4 scattering processes is negligible for any value of Tr, whereas the yield from 
the decay process q —* qa can dominate the abundance for Tr < rriq. 

The present fraction of the axino energy density to the critical density £l a is given by 

niaYa ~ 0.72 eV (8.1) 

where h is the Hubble parameter. Using this relationship, one can redisplay the results of 
Fig. in the (m a ,T R ) plane. 
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For axinos produced through scatterings and decays of particles in the plasma (TP 
axinos), eq. ( |8.1| ) gives an upper bound on Tr as a function of the axino mass from the 
requirement Qj p h 2 < 1. This bound is depicted by a thick solid line in Fig. ||, for m g = 



nig = 1 TeV and F = 10 11 GeV. 



At large Tr, the TP mechanism is dominant and the cosmologically favored region 
(0.2 < VL- a h 2 < 0.4) forms a narrow strip (not indicated in Fig. ||) just below the fiT p = 1 
boundary. 

For reheat temperatures less than the mass of the gluinos and squark, non-thermal 
production of axinos via out of equilibrium neutralino decays can dominate |13|, p^ ]. In 
that case the neutralinos decay into an axino and a photon or Z well after freeze-out, and 
the final axino abundance is independent of the reheat temperature, as long as this is high 
enough for the neutralinos to thermalize before freeze-out. The shaded region below the 
bound on Tr, shows where a sufficiently large axino population can be obtained by the 
decay of a nearly pure bino with mass m x = 100 GeV. The ranges of rria where axinos are 
hot/warm/cold dark matter are also shown. More details can be found in [12|. 

In Fig. ||, an enlarged region of the (ma,?^) plane is shown for different values of 
the sparticle masses. The dash-dotted line always represents the bound if squark decay is 
neglected; clearly, for heavy gluinos, a large region of the parameter space is excluded by 
including squark decay. This remarkably large effect is due to two factors: for one thing, 
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Figure 8: The (m^T^) plane, with a bino-like neutralino of mass 100 GeV, rrig — niq — 1 TeV 
and F = 10 11 GeV. The dash-dotted line represents the bound if squark decay is neglected, while 
the thick line is the full result. 

a light squark remains in thermal equilibrium and can play a role in axino production 
down to lower temperatures; for another, a large gluino mass enhances the squark decay 
to axinos. 

Note that, as discussed previously, the squark decay rate into axinos is only regulated 
by the ratio rrig/ f a . Therefore the same effect could be due to a smaller effective PQ scale 
f a rather than a larger gluino mass. 

9. Conclusions 

We have computed the effective axino-quark-squark vertex arising in KSVZ axion models 
through one loop diagrams, and we considered the effects of this interaction for the thermal 
production of axinos. 

As expected by dimensional arguments, we find that the dimension-4 scattering dia- 
grams generated by this coupling are negligible with respect to the dimension-5 scattering 
processes already analysed in On the other hand, the two body decay of a squark to 
an axino and a quark can be the dominant axino production channel at low reheat temper- 
ature, in particular when the squarks are substantially lighter than the gluinos. Our results 
put tighter constraints on the reheat temperature for the axino CDM scenario, summarized 
in Fig. ^. For very light squark masses and heavy gluinos, most of the parameter space for 
NTP axino DM is excluded. 
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Figure 9: The (ma,T#) plane, showing the different bounds on the reheat temperature depending 
on the mass spectrum. For light gluinos, the effective coupling is small and its effect on the axino 
abundance negligible. For light squarks, axino production through squark decay is so large that 
most of the allowed parameter space is excluded. Note that, in the figure in the left top panel, the 
kink at Tr ~ 10 2 GeV is where gluino decays start to dominate scatterings. For large gluino masses 
the kink is moved to the left and is less pronounced. In the right bottom panel, a similar kink is 
due instead to the decay of light squarks. 
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